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(54) FAST INTEGRATION METHOD AND SYSTEM 

(57)Abstract: 

PROBLEM TO BE SOLVED: To provide 
integration system/method effective even 
for highly dimensional Monte Carlo 
calculation, and high in convergence speed 
without impairing calculation time. 
SOLUTION: One-dimensional quantity is 
calculated from multidimensional random 
numbers. It is desirable that one- 
dimensional quantity is adjusted to the 
feature of a function to be integrated. 
One-dimensional quantity is classified into 
several groups by a value and the number 
of generation times is recorded at every 
group. The accumulation value of the 
value of the function to be integrated, which corresponds to the 
multidimensional random number at that time, is recorded in accordance 
with previous grouping. The operation is repeated for N times and a 
ratio between the distribution of one-dimensional quantity, which is 
logically predicted for the respective groups, and the recorded number 
of generation times is recorded as weight for the group. The product of 
the accumulation value of the integrated function and weight is 
calculated at every group and the sum of all the groups is taken. The 
sum is divided by the sum of weight against all the groups. Then, an 
integration value by corrected Monte Carlo analysis is given. 
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ro = 0.021673, 'a = 0.1817303, b = 0.0825398957, a 

void gettheory (double theory [STATNUM] ) ; 

double getaverage (double theory [STATNUM] , 

int us tat [STATNUM], double psum [STATNUM] , long cnt) ; 

int 

mainO 
{ 
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long 
double 
double 
int 



jj, maxcount = 10000; 

u[D]; 

norm, price, priceave; 

i, nindex, us tat [STATNUM] ; 
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double psumCSTATNUM], theory [ST ATNUM] ; 

for (i = 0; i < STATNUM; i++) 

{ 

ustat[i] = 0; 
psumCi] = 0; 

} 

.gettheory (theory) ; 
for (jj =1; jj <= maxcount;- jj++) 
■ { ■ 

getrandom(u) ; /* get D-ditn random sequences */ 

norm = calcnorm(u) ; /* calc norm of random sequence */ 
- nindex = (int) ( (norm+STATNUM*STATLEN / 2) / STATLEN) : 
price = value (u); /* get a price from a senario */ 
ustat [nindex]++; 
psum[nindex] += price; 

) 

priceave = getaverage (theory, ustat, psum, maxcount); 
printf ("result : %. IGfVn", priceave); 

} . ^ 

void gettheory (double theory [STATNUM] ) 

{ ■ ' 

int mid = STATNUM / 2; 
int i, j; 

double dx = STATLEN / TNUM; 
double c, f, sum, x, total = 0; 

c - dx / sqrt(2 * D * 3. 1415926536) ; , 
for (i = 0; i < mid - l; i++) 
{ - 
sum = 0; 

for (j = 0; j < TNUM; j++) 
{ 

x = i * STATLEN + (j + 0. 5) * dx; 
f = exp(- X * x / (2 # D))^ * c; 
sum += f; 

} 

theory [mid + i] = theory [mid - i -l] = sum; 
total += sum; 

} 

theory [0] = theory [STATNUM -1] = 0. 5 - total; 

} 

double getaverage (double theory [STAT^fUM] , 

int ustat [STATNUM], double psum [STATNUM], long cnt) 

( 

double pricesum = 0.0, wsum = 0.0, weight; 
int i ; 

for (i = 0; i < STATNUM; i++), 
{ 

weight = theory[i] / ustat[i]; 
pricesum += psum[i] * weight; 
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wsum weight * ustat[i]; 

) 

return pricesum / wsumi 
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